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ABSTRACT 

c~| ' We show that the condensates of a non-abehan gauge group, unified with the standard model 

gauge groups, can parameterize the present day cosmological constant and play the role of 
Qh. quintessence. The models agree with SNla and recent CMB analysis. 

' These models have no free parameters. Even the initial energy density at the unification scale 

and at the condensation scale are fixed by the number of degrees of freedom of the gauge 
^ ■ group (i.e. by Nc,Nf). The values of Nc,Nf are determined by imposing gauge coupling 

unification and the number of models is quite limited. Using Affleck-Dine-Seiberg superpotential 
one obtains a scalar potential V = A^+"</)~"'. Models with 2 < n < 4.27 or equivalently 
2 X 10~^GeV < Ac < 6 X W'^GeV do not satisfy the unification constrain. In fact, there are 
only three models and they have an inverse power potential with 6/11 < n < 2/3. Imposing 
primordial nucleosynthesis bounds the preferred model has Nc = 'i,Nf = 6, with n = 2/3, a 
condensation scale Ac = 4.2 x lO^^GeV and Wipo = —0.90 with an average value Wef f = —0.93. 
Notice that the tracker solution is not a good approximation since it has wtr = — = —0.75 
for n = 2/3. 

We study the evolution of all fields from the unification scale and we calculate the relevant cos- 
mological quantities. We also discuss the supersymmetry breaking mechanism which is relevant 
for these models. 
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1 INTRODUCTION 

The Maxima and Boomerang observations on the cosmic microwave background radiation 
("CMBR") and the superonovae project SNla have lead to conclude that the universe is flat 
and it is expanding with an accelerating velocity. These conclusions show that the universe is 
now dominated by a energy density with negative pressure with = 0.7 it 0.1 and tf^ < —2/3 
IQ]. New analysis on the CMBR peaks constrain the models to have w^o = — 0.82^'J^ 
This energy is generically called the cosmological constant. Structure formation also favors 
a non- vanishing cosmological constant consistent with SNla and CMBR observations An 
interesting parameterization of this energy density is in terms of a scalar field with gravitationally 
interaction only called quintessence S. The evolution of scalar field has been widely studied 



and some general approach con be found in p4|, 15]. The evolution of the scalar field (p depends 



on the functional form of its potential V{(j)) and a late time accelerating universe constrains the 



form of the potential |15|. 



It is well known that the gauge coupling constant of a non-abelian asymptotically free gauge 
group increases with decreasing energy and the free elementary fields will eventually condense 
due to the strong interaction, e.g. mesons and baryons in QCD. The scale where the coupling 
constant becomes strong is called the condensation scale Ac and below it there are no more 
free elementary fields. These condensates, e.g. "mesons", develop a non trivial potential which 
can be calculated using Affleck's potential |2^]. The potential is of the form V = A^^" 
where (p represents the "mesons", and depending on the value of n the potential V may lead 
to an acceptable phenomenology. The final value of u;^o (from now on the subscript "o" refers 



to present day quantities) depends n and the initial condition Jl^j [25]. A W(j)o < —2/3, which 
is the upper limit of |5[, requires n < 2.74 for Jl^j > 0.25 ]|2^. For smaller one obtains a 
larger tf^o for a fixe n. The position of the third CMBR peak favors models with n < 1 and 
for some class of models with V = M'^~^"'(f)~"'e'^'^^'^ , with n > l,/3 > 0, the constraint an the 
equation of state is even stricter — 1 < w^o ^ —0.93 ]0. In this kind of inverse power potential 
models (i.e. n < 2) the tracker solution is not a good approximation to the numerical solution 
because the scalar field has not reached its tracker value by present day. 

Here we focus on a non-abelian asymptotically free gauge group whose gauge coupling constant 



is unified with the couplings of the standard model ("SM") ones [13, 25]. We will call this 
group the quintessence or Q group. The cosmological picture in this case is very pleasing. 
We assume gauge coupling unification at the unification scale Ag^t for all gauge groups (as 
predicted by string theory) and then let all fields evolve. At the beginning all fields, SM and 
Q model, are massless and red shift as radiation until we reach the condensation scale Ac of 
Q. Below this scale the fields of the quintessence gauge group will dynamically condense and 
we use Affleck's potential to study its cosmological evolution. The energy density of the Q 
group 0,^ drops quickly, independently of its initial conditions, and it is close to zero for a long 
period of time, which includes nucleosynthesis (NS) if Ac is larger than the NS energy A^s (or 
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temperature T^s = 0.1 — WMeV), and becomes relevant only until very recently. On the other 
hand, if Ac < Atvs" than the NS bounds on relativistic degrees of freedom must be imposed on 
the models. Finally, the energy density of Q grows and it dominates at present time the total 
energy density with the — 0.7 and a negative pressure w^po < —2/3 leading to an accelerating 
universe 

The initial conditions at the unification scale and at the condensation scale are fixed by the 
number of degrees of freedom of the models given in terms of Nc, Nf. Imposing gauge coupling 
unification fixes Nc, Nf and we do not have any free parameters in the models (but for the susy 
breaking mechanism which we will comment in section ||). It is surprising that such a simple 
model works fine. 

The restriction on N^ Nf by gauge unification rules out models with a condensation energy 
scale between 2 x < Ac < 6 x lO^GeV or for models with 2 < n < 4.27 (the scale Ac is 



given in terms of Hq and n by Ac ~ ffo^^^^"^ [|ll|] , ^5|] ) . Since Wipo < — 2/3 requires n < 2.74 all 



models must then have Ac < 2 x 10~^GeV. The number of models that satisfy gauge coupling 
unification with a t(;^o < —2/3 is quite limited and in fact there are only three different models 
p5| . All acceptable models have n < 2/3 which implies that the condensation scale is smaller 
than the NS scale. The preferred model has A^'c = 3, Nf = 6, n = 2/3 and it gives w,f,o = —0.90 
with an average value Weff = —0.93 agreeing with recent CMBR analysis ^ |^. 

It is worth mentioning that we have taken Ac as the one loop renormalization energy scale (as 
used by Affleck et al l2l[] ) and if we had used the all loop renormalization energy scale the 
values of A'^c, Nf of the models may differ slightly but the general picture remains the same, i.e. 
there are only a few models that satisfy the requirement of gauge coupling unification, non of 
them have n > 2 and there are no free parameters. 



2 Condensation Scale and Scalar Potential 

We start be assuming that the universe has a matter content of the supersymmetric gauge groups 
SU{1) X SU{2) X SU{3) X SU{Q) where the first three are the SM gauge groups while the last 
one corresponds to the "quintessence group" Q and that the couplings are unified at Kgut with 
91 = 92 = 93 = 9Q= 9gut- 

The condensation scale Ac of a gauge group SU{Nc) with Nf (chiral + antichiral) matter fields 
has in = 1 susy a one-loop renormalization group equation given by 

Ac = Agute ^ (1) 

where bo = 3Nc — Nf is the one-loop beta function and Agut,ggut are the unification energy 
scale and coupling constant, respectively. From gauge coupling unification we know that Agut — 
10^6 GeV and ggut ^ V4vr/25.7 [||]. 
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A phase transition takes place at the condensation scale Ac, since the elementary fields are free 
fields above Ac and condense at Ac. In order to study the cosmological evolution of these con- 
densates, which we will call (j), we use Affleck's potential [21|. This potential is non-perturbative 
and exact |27|. 

The superpotential for a non-abelian SU{Nc) gauge group with Nf (chiral + antichiral) massless 
matter fields is [^] 

W={N,- Nf){ )yiN.-Nf) (2) 

det <QQ > 

where bo is the one-loop beta function coefficient. Taking det < QQ >= H^lifp] one has 
W = {Nc-Nf){Al°(j)~'^^f)'^/'^^'-'^f\ The scalar potential in global supersymmetry is y = jW^P, 



with = dW/dcj), giving |T], |T| 

V = c2A^+"0-" (3) 

with c = 2Nf, n = 2 + 4 ^y'^^^y and Ac is the condensation scale of the gauge group SU{Nc). 
The natural initial value for the condensate is (pi = Ac since it is precisely Ac the relevant scale 
of the physical process of the field binding. 

In eq.(^) we have taken cj) canonically normalized, however the full Kahler potential K is not 
known and for i;A ~ 1 other terms may become relevant |11| and could spoil the runaway and 
quintessence behavior of (p. Expanding the Kahler potential as a series power K = \(p\'^ + 
T,iai\4>\'^'' /2i the canonically normalized field 4>' can be approximated^ by (j)' = {K'^Y^'^cj) and 
eq.(|) would be given hy V = (xf = (2iVj)2A4+"0-"(i^J)("/2-i). For n < 2 the 



exponent term of is negative so it would not spoil the runaway behavior of (j) |^, 25]. 

If we wish to study models with < n < 2, which are cosmologically favored |25| we need to 
consider the possibility that not all Nf condensates (pi become dynamical but only a fraction u 
are (with Nf > v > 1) and we also need Nf > Nc |13, 25]. It is important to point out that 
even though it has been argued that for Nf > A'c there is no non-perturbative superpotential W 
generated ]21], because the determinant of QQ in eq.@ vanishes, this is not necessarily the case 
]23]. If we consider the elementary quarks Qf,Qf = 1,2, A^'j, a = 1,2, ...,Nc) to be the 
relevant degrees of freedom, then for A'^c < A^/ the quantity det{Q^oiQ'j) vanishes since, being the 
sum of dyadics, always has zero eigenvalues. However, we are interested in studying the effective 
action for the "meson" fields (/>*• =< Q^Q" >, and the determinant of i.e. det < Q^Q'j >, 
being the product of expectation values does not need to vanish when A^c < ]^f (the expectation 
of a product of operators is not equal to the product of the expectations of each operator). 

One can have u ^ Nf with a gauge group with unmatching number of chiral and anti-chiral fields 
or if some of the chiral fields are also charged under another gauge group. In this case we have 
c = 2v,n = 2 + 4 jy ^jy^ and Nf — v condensates fixed at their v.e.v. < QQ >= A^ Jl3]. Another 
possibility is by giving a mass term to A'^^ — v condensates =< QkQk >, (k = 1, •••,^/ ~ ^) 

^The canonically normalized field p' is defined as p' = g{p, p)p with Kt = {g + pg^ + pgsY 
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while leaving v condensates 0^ =< QjQj >, (j = 1, i^) massless. Notice that we have chosen 
a different parameterization for ip and (p. The mass dimension for 99 is 2 while for it is 1. The 
superpotential now reads, 

W = {N^- Nf){^J^f'^''^-''^^ + (4) 

with m the mass of QkQk- If we take the natural choice = A^, as discussed above, and m = 
|13| and we integrate out the condensates (p using 

dW 



dip 



we obtain ip = (Nj — zy)(^<:~^/)/{^c~'^)A^. By integrating out the ip field the second terms in 
eq.(^), which is proportional to the first term, can be eliminated. Substituting the solution of 
eq.(^) into eq.(^) one finds 

with a = 2v/{Nc - Nf). 

The scalar potential V = is now given by 

V = c'2A^+"'0-"' (7) 

with c'2 = 4i/2(^We^^2(^^_j,)(7V^-v)/(Afe-z.) ^ ^ 2 + 4zy/(iVc-A^/)- Notice that for v = ISI f 
we recover eq.@. From now on we will work with eq.(0) and we will drop the quotation on n! . 

The radiative corrections to the scalar potential eq.(0) are V ~ A^+"0~"(1 + 0(A^i;^~^)) 
They are not important because we have > Ac and are negligible at late times when ^ S> Ac. 



2.1 Gauge Unification Condition 

In order to have a model with gauge coupling unification the scale Ac given in eq.(^) or (^) 
must be identified with the energy scale in eq.(2). However, not all values of Ac will give 
an acceptable cosmology. The correct values of Ac depend on the cosmological evolution of the 
scalar condensate which is determined by the power n in eq. (0). The Ac scale can be expressed 
in terms of present day quantities by |^] 

Ac=(3i7o2yo2Cc-2)^ (8) 

where is the fraction of the total energy density carried in F, = y/3i/^ = ^^(l — w^jl^ 
and for = 0.7, w^o = —2/3 one has yo = 0.76. A rough estimate of eq.(P) gives Ac ~ 
j_j-2/{i+n) g-j^pg ^Yso expect (po = 0{1) today (we are living at the beginning of an accelerating 
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universe). The number of models that satisfy the unification and cosmological constrains of 
having fi^o = 0.7, /lo = 0.7 (with the Hubble constant given by Hq = lOO/io km/Mpc sec) and 
W(j,o < —2/3 1^] is quite limited [^5|. In fact there are only three models given in table 
These models are obtained by equating Ac from eq.(||), which is a function of Nc, Nf through 
bo, and eq.(|^, which is also a function of Nc, Nf, v through n. The exact value of yo, 't'o must 
be determined by the cosmological evolution of (p (c.f. eqs.(|l^)) starting at Ac until present day. 
For an acceptable model the parameters A'^c; A'^^ and v must take integer values. We consider 
an acceptable model when Ac in eqs.(||) and (^) do not differ by more than 50%. With this 



assumption there are only 3 models, given in table 2.1, that have (almost) integer values for 
Nf. In all these models one has n < 2/3 and the quantum corrections to the Kahler potential 
are, therefore, not dangerous. All other combinations of Nc,Nf,v do not lead to an acceptable 
cosmological model. 

From eq.(^) one has for n < 4.27 a scale Ac < 6.5 x lO'^GeV and from eq.(|^) this implies that 
bo < 5.7. Since bo = 3Nc — Nf = 2Nc + 4z^/(n — 2) and the minimum acceptable value for 
Nc is two one finds bo > 4: + Av/[n — 2). Taking 2 < n < 4.27 gives a value of bo > 5.7. 
The value of n = 4.27 gives the upper limiting value for which we can find a solution of eqs. 

and (^). We see that it is not possible to have quintessence models with gauge coupling 
unification with 2 < n < 4.27. In terms of the condensation scale the restriction for models with 
2 X IQ-^GeV < Ac < 6 X W^GeV. 

Using n = 2 + 4u/{Nc - Nf) or equivalently Nf = Nc + 4:u/{n - 2) with bo = 3Nc - Nf = 
2Nc + 4u/{n - 2) we can write from eq.(|l|) as bo = Sir^ /g'^g^t{Log{^))-^ and Nc 

Nc = -bo + 



2 2-n 

9gut ' 2-n 



-iLogi^-^])-' + ^ (9) 



Form eq.(|8|) we have Ac as a function of n (with the approximation of yo4>o — 1) in 
eq. (P) becomes a function of n and u only. In figure |l] we show A'^c as a function of n or Ac with 
the constraint of gauge coupling unification. We see that for 2 x W~^GeV < Ac < 6.5 x lO^Gey 
we have a A'c < 2 and therefore are ruled out. In terms of n the condition is that models with 
2 < n < 4.27 are not viable. In deriving these conditions, we have taken v = \ which gives the 
smallest constraint to A^c as seen from eq.(^). 

The upper limit Ac > 6.5 x lO^Ge^ has n > 4.27 (c.f. eq.(|8|)). As mentioned in the introduction. 



the value of w^o depends on the initial condition fl^i and on n . The larger n the larger u)^o 
will be (same is true for the tracker value wtr = —2/(2 + n)). It has been shown that assuming 
an initial value of no smaller than 0.25 then the value of ^i;^o will be less then < — 2/3 
only if n < 2.74 |^]. Therefore, the models with n > 4.27 are not phenomenological acceptable 
and since 4.27 > n > 2 are also ruled out by the constrain on gauge coupling unification, we are 
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Num 




Nf 


u 


n 


I 


3 


5.98 


1 


0.66 


II 


6 


14.97 


3 


0.66 


III 


7 


18.05 


4 


0.55 



Table 1: Models that satisfy gauge coupling unification and have n < 2.74 {i.e.w^o < —2/3) 
left with models with 

Ac < 2 X lO'^GeV or n < 2. (10) 
So, only models with a cosmological late time phase transition are allowed. 




2 4 6 8 10 -10 -5 5 10 

n Log^o^c GeV 



Figure 1: We show Nc as a function of n and the energy scale Ac after imposing gauge coupling 
unification. Nc must be larger than 2 and we have taken 1^ = 1. 



3 Thermodynamics, Nucleosynthesis Bounds and Ini- 
tial Conditions 

Before determining the evolution of (p we must analyze the initial conditions for the SU{Q) gauge 
group. The general picture is the following: The Q gauge group is by hypothesis, unified with 
the SM gauge groups at the unification energy A^^^. For energies scales between the unification 
and condensation scale, i.e. Ac < A < Agut, the elementary fields of SU{Q) are massless and 
weakly coupled and interact with the SM only gravitationally. The Q gauge interaction becomes 
strong at Ac and condense the elementary fields leading to the potential in eq.(0). 

Since for energies above Agut we have a single gauge group it is naturally to assume that all 
fields (SM and Q) are in thermal equilibrium. However, at temperatures T < Tg^t the gauge 
group Q is decoupled since it interacts with the SM only via gravity. 
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The energy density at the unification scale is given by pTot = ^9TotT^, where gxot = T,Bosons + 
7/8TiFermions is the total number of degrees of freedom at the temperature T. The minimal 
models have grot = Qsmi + QQi, with Qsmi = 228.75 and qqi = (1 + 7/8){2{N^ - 1) + '^NjNc) for 
the minimal supersymmetric standard model MSSM and for the SU{Q) supersymmetric gauge 
group with Nc colors and Nf (chiral + antichiral) massless fields, respectively. The initial energy 
density at the unification scale for each group is simply given in terms of number of degrees of 
freedom, Q = p/pc, 

i^Qi{Agut) = , -Bcm VLsmii^gut) = (H) 

gTot QTot 

with = + = 1- Since the SM and Q gauge groups are decoupled below Kgut, their 
respective entropy, Su = gua^T^ with the degrees of freedom of the k group and a the scale 
factor of the universe (see eq.(|20|)), will be independently conserved. The total energy density 
p as a function of the photon's temperature T above Ac (i.e. Ac < A < Agut), with the Q fields 
still massless and redshifting as radiation, is given by 



P=^..T' (12) 



with 

4/3 



r 



g* = gsrnf + gQf(^] = gsmf + ggf [ ^f^^ 1 (13) 
V J / \gsmigQf J 

and gsmi,gsmf,gQi,gQf are the initial (i.e. at the unification scale) and final standard model 
and Q model relativistic degrees of freedom, respectively. From the entropy conservation, we 
know that the relative temperature between the standard model and the Q model is given by 
T'/T = (££HL££2i)i/3, It is clear that the energy density for the Q model pg = Tr^/SOgqT'^ in 
terms of the photon's temperature T is fixed by the number of degrees of freedom. 



n 



Qf 



g*T^ 

srn fgQi/gsmigQfY^'^ 
gsmf + gQf{gsmfgQi/ gsmigqfY/^ 



(14) 



Eq.([lJ) permits us to determine the energy density of the Q group at any temperature above 
the condensation scale. 



3.1 Energy Density at 

We would like now to determine the energy density at the condensation scale which will set the 
initial energy density for the scalar composite field (p. 

Just above the condensation scale Ac we take, for simplicity of argument, that all particles in 
the Q group are still massless and we can use eq.(14) to determine the fig (Ac) with gqi = ggf 
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giving 

gQfiSsmf / dsmi) _n 

= ] 7 -/ W3- U^) 

dsmf + 9Qf[9smf/9smi) ' 

At Ac we have a phase transition and we no longer have elementary free particles in the Q group. 
They are bind together through the strong gauge interaction and the acquire a non-perturbative 
potential and mass given by eq.(0). In other words, below the condensation scale there are no 
free "quarks" Q and we have "meson" and "baryon" fields. 

If we consider only the SM and the Q group, the energy density within the particles of the Q 
group must be conserved since they are decoupled from the SM (the interaction is by hypothesis 
only gravitational). 

Furthermore, the "baryons", which we expect to be heavier then the lightest meson field (as in 
QCD), and the massive "meson" fields ip (see eq.(Q)) are coupled (i.e. T/H > 1) to the lightest 
composite field (p for temperatures Agut > T > Ac(Ac/mp)^/^, with Ac/rup <^ 1. The "baryons" 
and the heavy ip fields will then decay into the lightest state within the Q group, i.e. the (p field. 
So, we conclude that all the energy of the Q group is transmitted into (p at around the phase 
transition scale given by condensation scale Ac and 

n^Ac) = flQiAc). (16) 

This is a natural assumption from a particle point of view but is not crucial from a cosmological 
point of view, in the sense that any " reasonable" fraction of the energy density in the Q group 
would give a correct cosmological evolution of the (j) field. 

We would like to stress out that the initial condition for (f> is no longer a free parameter but it 
is given in terms of the degrees of freedom of the MSSM and the Q group. 



3.2 Nucleosyntehsis Constrain on 

The big-bang nucleosynthesis (NS) bound on the energy density from non SM fields, relativistic 



or non-relativistic, is quite stringent < 0.1 [19| and a recent more conservative bound gives 
no < 0.2 i2C 



If the Q gauge group condense at temperatures much higher than NS then, the evolution of the 
condensates will be given by eqs.(^) with the potential of eq.(^) and we must check that Jig 
at NS is no larger than 0.1-0.2. This will be, in general, no problem since it was shown that 
even for a large initial Og at the condensation scale the evolution of (j) is such that fig decreases 
quite rapidly and remains small for a long period of time (see figure |2|) [^, 25]. 



On the other hand, if the gauge group condenses after NS we must determine if the Q energy 
density is smaller than J7g < 0.1 — 0.2 at NS. Since the condensations scale Ac is smaller than 
the NS scale, all fields in the Q group are still massless and the energy density is given in terms 
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of the relativistic degrees of freedom and from eq.(|T^ to set a limit on qqj and ggi, 

A _ -1/3 4/3 -1/3 4/3 

^9Q = 9Qf 9Qi = 9smf 9smi (17) 

and for qqj = gqi = gq 

A -1/3 4/3 

^9Q =9Q = Y~Q~3smf 9smi (18) 

where we should take gsmf = 10.75 at the final stage (i.e. NS scale) and gsmi = 228.75 at 
the initial stage (i.e. at unification) for the minimal supersymmetric standard model MSSM. 
For fig < 0.1,0.2 eq.(|T7|) gives un upper limit on the number of relativistic degrees of freedom 
Ac/Q < 70, 158 respectively (or gq < 70, 158 if gqj = gqi = gq). 

The l.h.s. of eq.(|l7|) depends on the initial (i.e. at unification) and final (at NS) number of 
degrees of freedom of the gauge group Q. The smaller (larger) the initial (final) degrees of 
freedom of Q the smaller Agq and 0,q will be. 



3.3 Supersymmetry Breaking 

Another important ingredient in these models is the way supersymmetry is broken. The precise 
mechanism for susy breaking is still an open issue but it is generally believed that gaugino 
condensation of a non-abelian gauge group breaks susy There are a two ways that the 



breaking of susy is transmitted to the MSSM, by gravity [^] or via gauge interaction |3C]. 

In the case of gravity susy breaking, the same mechanism that breaks susy for the MSSM will 
break susy for the Q group and from particle physics we expect the breaking to be transmitted 
at m ~ 

Hreak/'^l ~ TeV scale (i.e. Abreak ^ W^^GeV). The final de grees of freedom of 
the Q group must contain only the non-super symmetric ones at temperatures T < TeV, with 
gqj = 2{N^-l) + 2NfNc7/8 at NS and the initial ones at unification are gqi = {1 + 7 /8){2{N^ - 
1) + 2NfNc). The Q group would be globally supersymmetric but would have explicit soft 
supersymmetry breaking terms (as the breaking of MSSM to SM). The fields in the gauge group 
responsible for susy breaking are not in thermal equilibrium at T < Tg^t neither with the SM 
nor with the Q group since they interact via gravity only. 

On the other hand if susy breaking is gauge mediated and since the Q group interacts only 
gravitationally with all other gauge groups, the supersymmetry breaking for Q group will be at a 
scale m ~ A^ireak/^p ~ W~^^GeV , since one expects the condensation scale of the susy breaking 
gauge group to be in this case much smaller than for the gravity one, with Abreak < 0{W'^GeV) 
|l30|| , to give a susy breaking mass to the SM of the order of TeV. Therefore, in this second case 
the Q group will be supersymmetric for models with Ac > m ~ lO^^^Ge^ and the relativistic 
degrees of freedom at NS will be the same as the initial ones, i.e. gqf = gqi = (1 + 7/8){2{N^ — 
1) + 2NfNc) at NS. If susy breaking is gauge mediated, then the gauge group responsible for 
susy breaking will be coupled to the MSSM and will be in thermal equilibrium at Ac < T < Tgut 
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and its degrees of freedom must be taken into account in the initial gsmi = gsAU + 9ex, where 
9s Mi are the degrees of freedom of the MSSM and Qex those of the gauge group responsible 
for susy breaking. Typical models of susy breaking via gauge interaction have a gauge group 
SU{Nc) with Nc> 5 and iV/ > 4 @ which gives and extra gex > 160. 

We would like to point out that in both cases the susy breaking mass is a problem for quintessence 
since the present day mass must be of the order of W~^^GeV, much smaller than the susy 
breaking mass. Here, we have nothing new to say about this problem and we consider it as 
part of the ultraviolet cosmological constant problem, i.e. the stability of the vacuum energy 
(quintessence energy) to all quantum corrections. The contribution to the scalar potential from 
the susy breaking scale from the 4> field and/or from any other field of the MSSM is enormous 
compared to the required present day value. The ultraviolet problem is an unsolved and probably 
one the most important problems in theoretical physics. 



3.4 Models 

Now, let us determine the contribution to the energy density at NS for the three models given 
in table |2|, taking the closest integer for Nj. The number of degrees of freedom for an SU{Nc) 
supersymmetric gauge group with Nf flavors is gqi = (1 + 7/8){2{N^ — 1) + 2NfNc). All 
three models have the same supersymmetric one-loop beta function bo = 3Nc — Nf = 3 and 
Ac = 4 X IQ-^GeV. 

The group with the smallest number of degrees of freedom is Model I, Nc = 3, A'^^ = 6 and we have 
dQs = 97.5 supersymmetric degrees of freedom. In this model, we have at NS 0,q\ns = 0.13. 
We see that the energy density of Q is slightly larger than the stringent NS bound ^q\ns < 0.1 
but it is ok with the more conservative bound 0,q < 0.2. 

For other thwo groups we have ^Iq\ns = 0.42,0.51 for Models II and III respectively. ^q\ns is 
larger since they have a larger gqi (i.e. Nc,Nf) and all these models would not satisfy the NS 
energy bound ^q\ns < 0.2. Therefore, if susy is broken via gravity these three models would 
not be phenomenologically viable, unless more structure is included. On the other hand, if susy 
is broken via gauge interaction we would need to take into account the degrees of freedom of the 
susy breaking group given in gsmi = dSMi + dex when calculating ^q\ns- These extra degrees 
of freedom give a larger gsmi and therefore reduce Q^i as can be seen from eq.(14). In order to 



have ^q\ns ^ 0.1 we require gex ^ 64, 718,986 for models I, II, and III, respectively, while for 
^q\ns ^ 0.2 we require gex ^ 287,433 for models II and III, respectively. 

We have checked that a larger number of extra degrees of freedom does not affect the cosmo- 
logical evolution of 4> significantly. In fact there is no "reasonable" upper limit on gex from the 
cosmological point of view (e.g. for gex = 10^ the model is still ok) as can be seen in fig.§. Notice 
that a large gex ^ 1 gives a small energy density Q^i oc gsmi^ oc Qex^^ ■ This result also shows 
that an acceptable cosmological model cosmological is almost independent on the initial energy 
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density of 4>. 

As a matter of completeness we give the minimal model when susy is broken via gravity. In this 
case one has to consider that Qqj QQi as discussed in section 3^. The minimal gauge group 
when susy is broken via gravity has Nc = 5, Nf = 14, i/ = 1 and qqs = 352.5, qqus = 170.5 
for the relativistic susy and non-susy degrees of freedom, respectively, and one has n = 14/9 ~ 
1.5, Ac = 4.5 X W~^GeV, ^q\ns = 0.41 much larger than the NS bound. The difference in the 
values of Nc, Nj between the susy and non-susy models are due to a change in bo, the one loop- 
beta function in eq. (I), below the susy breaking scale ITeV, giving different values for Ac for the 
same A^'c, Nf. We conclude that unless more structure is included (i.e. need gex = 689 relativistic 
fields coupled to the SM to have ^q\ns < 0.1) there are no models that satisfy the NS energy 
bound for the case when susy is broken via gravity. However, if we allow for a discrepancy in 
Ac from eqs.(|l|) and of up to one order of magnitude then the model Nc = 3, Nj = 9,v = 3 
would be fine and it has qqs = 131.25, gQns = 63.25 for the relativistic susy and non-susy 
degrees of freedom, respectively, with n = 4/3, Ac = 1 — 9 x lO^^^GeV, 0.q\ns = 0.2. 



4 Cosmological Evolution of 



The cosmological evolution of (j) with an arbitrary potential V{(f)) can be determined from a 
system of differential equations describing a spatially flat Friedmann-Robertson- Walker universe 
in the presence of a barotropic fluid energy density that can be either radiation or matter, 
are 

P = -3H{p + p), (19) 



where H is the Hubble parameter, cj) = dcp/dt, p (p) is the total energy density (pressure). We 
use the change of variables x = and y = and equations (|l9|) take the following form 
11,11]: 

-A + -x[2x^ + 7^(1 -x^- y^)] 

VN = -^Xxy + ^y[2x^ + -f^{l-x^-y^)] (20) 
Hn = -Ih[2x^+j^{1-x^ -y^)] 



where N is the logarithm of the scale factor a, N = ln{a); fj\f = df/dN for / = x,y,H; 
7^ = l + w^ and A(A^) = —V'/V with V = dV/dcp. In terms of x, y the energy density parameter 
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2 2 

is Q(f, = + while the equation of state parameter is given by Wtj, = = ^^iZ^ji (with 



ml = G/Svr = 1). 

The Friedmann or constraint equation for a flat universe fi-^ + = 1 must supplement equations 
( pO| ) which are valid for any scalar potential as long as the interaction between the scalar field 
and matter or radiation is gravitational only. This set of differential equations is non-linear 
and for most cases has no analytical solutions. A general analysis for arbitrary potentials is 
performed in |jl^, the conclusion there is that all model dependence falls on two quantities: 
A(A^) and the constant parameter 7-y. In the particular case given by ^ tx we find A — > 
in the asymptotic limit. If we think the scalar field appears well after Planck scale we have 
Aj = nmpi/cpi = nmpi/Ac ^ 1 (the subscript i corresponds to the initial value of a quantity). 
An interesting general property of these models is the presence of a many e-folds scaling period 
in which A is practically a constant and <C 1 . After a long permanence of this parameter at a 



constant value it evolves to zero, A — > 0, which implies ^ < and ^ > |15|, leaving us with 

= + ^ 1 and uJj^o = ^T^p' ~^ ~1) which are in accordance with a universe dominated 
by a quintessence field whose equation of state parameter agrees with positively accelerated 
expansion. 

The evolution of fi^ can be observed in Figure together with the evolution of which fulfills 
the condition < —2/3 [Q] for different initial conditions. 

The value of the condensation scale in terms of Hq is 

Ac = [-^j <^t-- (21) 

and together with eq.(|^) sets a constrain for Nc,Nf. The approximated value for y1,(f>o can 
be obtained from eq.(p3|) but one expects in general to have 0.76 < Vo < 0.83 and i?!io ~ 1 for 
= 0.7 and w^o < —2/3. This can be also seen from the identity = 0(^(1 — w^)/2. The 
order of magnitude of the condensation scale is therefore Ac = Ho^^^^""^ . 

The value of Wcj,o can be approximated by |25] 

«^0o = -l + ^r# 22 



with given by solving |25| 



l%o = ^ (23) 



where (/>sc is the scaling value of i.e. the constant value at which </> stays for a long period 
of time. The scaling value is given only in terms of O,^,, (^^c = ^Pi + ^/Q^XjA for < 1/2 and 
(psc = ^i + V6{j^ + iLo(?[^]) for % > 1/2 i. 

In order to analytically solve eqs.(^) we need to fix the value of n and we can determine 



w^o by putting the solution of (p3D into eq.(|24). Eq.(23) can be rewritten as 0o = </'sc(l — 
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'n?^(f,o/Q(t)'i) ^-^d we see that 4>o > <psc- For 7^0 = n'^^}^o/^<Pl ^ 1 one has (f)o — 4>sc and 
for the simple cases of n = 1,2 and 4 we find (/>o|ra=i = + \J^4''ic + GJl^o/S, yoU=i = 

<^sc(-3./..c + y9^2~^61V> ) <^o|n=2 = y^2~P2?v73, y2|^^2 = 3</.2^1^^„/(3(^2^ + 20<^o) and 

0o|n=4 = ^40,^0/3 + y^902^ + lerj^^o/S, y2|^^4 = - 8J^^„/(4f^^o + ^g^^TiefV), respec- 
tively. Notice that the value of (pojW^po at fi^o = 0.7 cioes noi depend on Hi or iJo and it only 
depends on (through (psc) and n. 



5 The Models 
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Figure 2: We show the evolution for ^^^w^ for initial condition = 0.07 dashed-dotted and 
dashed lines, respectively and for Vl^i = 10~^°, dotted and solid lines, respectively, with n = 2/3. 
The first case corresponds to gex = 64 while the later case has a huge number of extra degrees 
of freedom g^^ > 10^. The vertical lines correspond to present day values with = 0.7 and 
K = 0.7. 



In this section we study the three different models given in table |2.1| . It is interesting to note 
that all three models have a one-loop beta function coefficient bo = 3Nc — Nj = 3 which implies 
that they have the same condensations scale Ac = 4.2 x 10~^GeF. The power of the exponent 



n, see table is very similar and if we take the closest integer value for A'^^ one has n = 2/3 
or 6/11. Notice that model I is self dual Nc = Nf — Nc = 3 with Nj matter fields. The other 
two models are not self dual. 

From now on we will focus on the Model I of table § and we will summarize the relevant quantities 
in table H for all models. 



The initial energy density at the unification scale is given by eq.(14) with gqi = 97.5, gsmi = 
228.75 is Qq^Agut) = 0.3. Below Agut the fields are weakly coupled, massless (they redshift 
as radiation) and are decoupled from the SM. A phase transition takes place when the gauge 
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coupling constant becomes strong at the condensation scale. Since the condensation scale is much 
smaller than the NS scale, Ac <C O.lMeV, we expect all fields of the Q group to be relativistic at 
NS. From eq.(|l|) with 

9smf — 10.75 the energy density, assuming no extra degrees of freedom, 
is ^q\ns = 0.13 for susy Model I. In order to satisfy the NS bound ^q\ns < 0.1, 64 extra 
relativistic degrees of freedom in thermal equilibrium with the SM at T < Tg^t are required 
while for r^glArs < 0.2 the model does not require any extra degrees of freedom. 

What is the energy density of Q at the condensation scale A^? Using eq. (|14| ) with (jsmf — 
3.36, 5q/ = 97.5 given at Ac = 4.2 x lO^^Gey and with no extra degrees of freedom at the 
unification scale (i.e. gsmi = 228.75 and gqi = gqj = 97.5) for Model I one has Q.^{A.c) = 
Oq(Ac) = 0.095. Imposing the stringent NS bound ^q\ns < 0.1 we need to include g^x — 

64 

extra degrees of freedom (that should come from susy breaking mechanism) and the energy at 
Ac is now Vt^{Kc) = 0.07. 

Evolving eqs.(]20|) with initial condition 0(^(Ac) = 0.07 gives at present time with ho = 0.7, 0,^0 = 
0.7 a value of w^ff = J da Q{a)w^{a)/ J da Q{a) = —0.93 with of w^o = —0.90 in agreement 
with SNla and CMBR data. The analytic solution given in eq.(p^) is W(j,o{Th) = —0.82 and it is 
a much better approximation to the numerical value then the tracker value wtr = —2/(2 + n) = 
-0.75 [| which is the upper value of tf^o for given n and arbitrary initial conditions. 

From a cosmological evolution point of view, we have a large range of initial condition of ^^pi [ p^ . 
The upper limit is set by NS and there is no "reasonable" lower limit (a smaller implies that 
we have a much larger number of extra degrees of freedom g^x but it must be finite) still gives 
an acceptable model and there is clearly no fine tuning in these models. The effect of a large 
number of extra degrees of freedom gex ~ 10'^ at the condensation scale is to drop the energy 
density from Q(j,{Ac) = 0.07 with g^x = 64 to r2^(Ac) = 0.01 with gex = 10^ and the numerical 
solution, in this case, gives w^ffo = —0.82 at present time still within the observational limits. 
In fact, there is no upper limit for g^x from the cosmological evolution constrains for (p because 
the upper value for w^o is given by its tracker value which for n = 2/3 is wtr = —0.75 < —2/3 
smaller than the upper limit given by SNla and CMBR data. In figure |2| we show the evolution of 
Cl^, Wfj, for the minimal number of gex = 64 (Jl^j = 0.07) and for an extreme case with gex ~ 109 
(r^^j = 10^^*^) and in both cases we get an acceptable model. 

In table |2| and ^ we summarize the relevant cosmological quantities. In table |2| we give the 
values of n,bo, the deg rees of freedom of Q with {gQs) and without supersymmetry (^gQnsusy) ^ 
the condensation scale Ac. Notice that all models have same 6o, Ac but n differs slightly. Model 
I is the minimal model, in the sense that it has the smallest number of degrees of freedom. 

In table ^ we give the values of the initial energy density Q.Q{Kgut)t the energy density at NS 
(for gex = 0), the number of extra degrees of freedom needed to have Qq{NS) = 0.1 or 0.2, 
the value of Q^{Ac) with gex = 0, the value of Nxot (the e- folds from Ac to present day), the 
values of w^po and Weff calculated numerically and the value obtain analytically from eq.(|2^) 
gives a good approximation to the numerical one. The energy density at Ac with the condition 
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Num 


Nc 




u 


n 






9qs 


I 


3 


6 


1 


2/3 


3 


4.2 X 10~^ 


97.5 


II 


6 


15 


3 


2/3 


3 


4.2 X 10-8 


468.5 


III 


7 


18 


4 


6/11 


3 


4.2 X 10-8 


652.5 



Table 2: We show the matter content for the three different models and we give the number 
of degrees of freedom for the susy and non susy Q group in the last two columns, respectively. 
Notice that the condensation scale and ho is the same for all models. 



Num 




^q{NS) 




Qexl 


gex2 








NtoI 


I 


0.30 


0.13 


0.09 


64 





-0.90 


-0.93 


-0.82 


12.9 


II 


0.67 


0.42 


0.33 


718 


287 


-0.90 


-0.93 


-0.82 


12.9 


III 


0.74 


0.50 


0.41 


986 


433 


-0.93 


-0.95 


-0.87 


11.1 



Table 3: The first column gives the model number. In columns 2-4 we give the energy density at 
different scales assuming no extra degrees of freedom (i.e. gex = 0). In column 5 and 6 we show 
the necessary number of gex to have Q.q{NS) < 0.1,0.2, respectively. We show in columns 7, 8 
the present day value of Wfj, calculated numerically and in column 9 the theoretically obtained 
from eq.(^). Finally, we give in the last column the number of e- folds of expansion from Ac to 
present day. 

ngiNS) < 0.1 gives = 0.07 for all three models while for nq^NS) < 0.2 gives = 0.15 
for models II and III, respectively. Since the number of gex for models II and III is quite large 
(larger than MSSM) we consider them less "natural" then the minimal Model I. 

6 Summary and Conclusions 

We have shown that an unification scheme, where all coupling constants are unified, as predicted 
by string theory, leads to an acceptable cosmological constant parameterized in terms of the 
condensates of a non-abelian gauge group. These fields play the role of quintessence. 

Above the unification scale we have all fields in thermal equilibrium and the number of degrees 
of freedom for the SM and Q model determines the initial conditions for each group. Below Agut 
the Q group decouples, since it interacts with the SM only through gravity. For temperatures 
above the condensation scale of the Q group its fields are relativistic and red shift as radiation. 
The entropy of each systems is independently conserved and we can therefore determine the 
energy density at NS and at Ac. The models we have obtain have a condensation scale below 
NS and in order not to spoil the NS predictions the energy density must be ^Iq{NS) < 0.1 — 0.2. 
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Without considering the contribution from the susy breaking sector, all models have Qq > 0.1 
with the smallest contribution from Model I {N^ = 3, Nf = 6 minimal model) giving ^Iq = 0.13 
in agreement with the conservative bound ^q\ns < 0.2 NS bound and slightly large then 
^q\ns < 0.1. If susy is transmitted via gravity we require extra structure to agree with the 
strongest NS (the gauge group responsible for susy breaking is not in thermal equilibrium with 
the SM below Ag^t) but if susy is gauge mediated than the NS bound is alleviated since one has 
extra degrees of freedom in thermal equilibrium with the SM. The cosmological evolution of 
quintessence is not sensitive to the number of the extra degrees of freedom. There is a minimum 
number required from NS bounds but there is no upper limit. 

At the condensation scale the Q fields are no longer free and they condense. We use Affleck's 
potential to parameterize the condensates and we study the cosmological evolution with the 
initial condition determined in terms of Nc, Nf only. Gauge unification determines the values of 
Nc, Nf and there are no models with 2 x < Ac < 6 x lO^GeV or 2 < n < 4.27. Since 

w^o < —2/3 requires n < 2.74 all models must have Ac > 2 x lO^^GeV or n > 2. The three 
acceptable models have a potential of the form V ~ i?i>~" with 6/11 < n < 2/3. The value of n 
and the energy density at Ac determines the present day value of W(f,o- We show that the models 
have 0<^o = 0.7, = —0.90 with a Hubble parameter ho = 0.7 and the value of n and w^o 
are in accordance with constrains from recent CMBR analysis, i.e. n < 1 and it;<^o = — 0.82^'^'J^ 
1^, 1^ . We also show that the tracker solution to inverse power potential is not specific enough 
(in Model I wtr = —0.75) and does not give a good approximation for models with n < 2, which 
are the cosmologically favored. 

We would like to stress out that there are no free parameters, not even the Q initial energy 
density at unification nor at the condensation scale. The models are well motivated from a 
particle physics point of view, they involve a late time phase transition, and they agree with 
present day observations. 

This work was supported in part by CONACYT project 32415-E and DGAPA, UNAM project 
IN-110200. 
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